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Abstract 

We have investigated the isoscalar giant resonances in the Sn isotopes using inelastic scattering 
of 386- MeV a-particles at extremely forward angles, including 0°. We have obtained completely 
"background-free" inelastic-scattering spectra for the Sn isotopes over the angular range 0°-9° 
and up to an excitation energy of 31.5 MeV. The strength distributions for various multipoles were 
extracted by a multipole decomposition analysis based on the expected angular distributions of the 
respective multipoles. We find that the centroid energies of the isoscalar giant monopole resonance 
(ISGMR) in the Sn isotopes are significantly lower than the theoretical predictions. In addition, 
based on the ISGMR results, a value of K T = —550 ± 100 MeV is obtained for the asymmetry 
term in the nuclear incompressibility. Constraints on interactions employed in nuclear structure 
calculations are discussed on the basis of the experimentally-obtained values for and K T . 
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I. INTRODUCTION 



Isoscalar giant resonances have been extensively studied since the discovery of the 
isoscalar giant quadrupole resonance (ISGQR) in the early 1970s jl s| . The isoscalar giant 
monopole resonance (ISGMR)was identified in 1977 0, [5] and was the subject of a number 
of studies through the 1980s 6r|8j. The isoscalar giant dipole resonance (ISGDR) was first 
reported by Morsch et al. [9] in 208 Pb but was conclusively identified by Davis et al. 10 ]. 
Both ISGMR and ISGDR are classified as compression modes and provide information about 
nuclear incompressibility, Ka, from which the incompressibility of infinite nuclear matter, 
Kqo, may be obtained [11 ]. 

Most of the earlier investigations of the isoscalar giant resonances used inelastic a scat- 
tering at 100-200 MeV and the strength of a particular giant resonance was assumed to be 
concentrated in a single peak with a Gaussian or Lorentzian shape. The resonance parame- 
ters were obtained by multiple-peak fits to the inelastic scattering spectra, after subtraction 
of a suitable "background" |l2|, 



Q. In the last decade, the Texas A&M (TAMU) group 
has carried out (a, a') studies of many nuclei at a bombarding energy of 240 MeV and ex- 
tracted the strength distributions of various isoscalar giant resonances in a number of nuclei 
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24] using a multipole decomposition analysis (MDA) [25[. Contemporaneously, we have 



carried out giant resonance measurements u sing inelastic scattering of 386 MeV a particles 
at extremely small angles, including 0-0. An especially useful feature of our mea- 
surements has been the elimination of all instrumental background events from the inelastic 
scattering spectra which was rendered possible by the optical properties of the Grand Raiden 



spectrometer 



36|. 



Here, we report on measurements of the isoscalar giant resonances in the even-A Sn 
isotopes (A=112-124) over the excitation-energy range 8.5-31.5 MeV. Previously, giant res- 
onance measurements on the Sn isotopes have been reported by the TAMU group [?], 
and KVI group 8J, using inelastic a scattering at 120-130 MeV and peak- fitting analyses of 
spectra. More recently, the strength distributions of various isoscalar resonances have been 
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2l|. 



obtained in some Sn isotopes (A=112,116,124) by the TAMU group 

High-quality measurements of the ISGMR over the full range of isotopes provide the 
opportunity to investigate the asymmetry term, K T , of the nuclear incompressibility. This 
term, associated with the neutron excess (N — Z), is crucial in obtaining the radii of neutron 
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stars in equation-of-state (EOS) calculations [38l44ll]; the asymmetry ratio, [(N — Z)/A] 



changes by more than 80% over the range of the investigated Sn isotopes. 

II. EXPERIMENTAL TECHNIQUES 

The experiment was performed at the ring cyclotron facility of the Research Center 
for Nuclear Physics (RCNP), Osaka University, using inelastic scattering of 386-MeV a 
particles at extremely forward angles, including 0°. A 4 He ++ beam was accelerated by 
the Azimuthally Varying Field (AVF) cyclotron up to 86 MeV, injected into the K = 400 
ring cyclotron for acceleration up to 386 MeV, and achromatically transported to the WS 
experimental hall without any defining slits. To reduce the background at and near 0°, the 
beam halo has to be tuned carefully in the experiment. This was accomplished by tuning 
the beam profile of the injection beam for the ring cyclotron, and typically less than 1 out 
of 10000 events had contamination from the other bunches. The beam current was 1-20 nA, 
which was limited by the data acquisition rate and by the maximum available current of 
the accelerator. The energy resolution obtained was ~ 150 keV full width at half maximum 
(FWHM). 

Self supporting target foils of enriched even- A 112 " 124 Sn isotopes of thickness 5.0-9.25 
mg/cm 2 were employed; we used special target frames with a large aperture in order to reduce 
the background caused by the beam-halo hitting the frames. Data were also taken with a 
nat C target at the actual field settings used in the experiments and the energy calibration 
was obtained from the peak positions of the 7.652- and 9.641-MeV states in the 12 C(a, a') 
spectra. 

Inelastically-scattered a particles were momentum-analyzed with the high-resolution 



magnetic spectrometer, "Grand Raiden" 36j, and the vertical and horizontal positions of the 



a particles were measured with the focal-plane detector system comprised of two position 



sensitive multi-wire drift chambers (MWDCs) and two scintillators 29]. The MWDCs and 
scintillators enabled us to make particle identification and to reconstruct the trajectories of 
the scattered particles. The scattering angle at the target and the momentum of the scat- 
tered particles were determined by the ray-tracing method. The vertical-position spectrum 
obtained in the double- focusing; mode of the spectrometer was exploited to eliminate the 



instrumental background 28|, |29] 
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Giant-resonance data were taken with the spectrometer central angle (6 spec ) set at 0°, 
2.5°, 3.5°, 5.0°, 6.5°, and 8.0°, covering the angular range from 0° to 9.0° in the laboratory 
system. The actual angular resolution of the MWDCs, including the nominal broadening 
of scattering angle due to the emittance of the 4 He ++ beam and the multiple Coulomb- 
scattering effects, was about 0.15° [42]. The vertical acceptance was limited to ±20 mr by a 
2-mm-thick tantalum collimator. The energy bite of the Grand Raiden spectrometer and the 
special MWDC arrangements for the 0° measurements limited the excitation energy range 
observed to £L=8.5-31.5 MeV. 



The incident 4 He ++ beam was stoppec 
to the different settings of Grand Raiden 



at three independent Faraday Cups (FC) according 



291 ] . In the measurements at large angles (6 > 6.5 C 



the beam was stopped on the FC mounted inside the scattering chamber (SC-FC). For 



291 ] ; the incident 4 He 



++ 



measurements at 0°, the FC was located downstream of the MWDCs 
particles were guided to a vacuum pipe situated at the high- momentum side of the MWDCs 
and finally stopped at the 0°-FC. A third FC was used for measurements in the scattering- 
angle region 2° < 9 < 5°. This FC was installed behind the Ql-magnet of the Grand 
Raiden (Ql-FC) 43[]. The use of these three Faraday-cups allowed us to obtain reliable 
values of accumulated charges for the incident 4 He ++ beam at different scattering angles. 
Normalization of the FCs was obtained with a beamline polarimeter located upstream from 
the target. The polarimeter target was inserted in the beam, and the scattered 4 He ++ 
particles counted, before and after each change of the FC. The overall accuracy of this 
normalization is estimated to be ~2%, including systematic errors from slight changes in 
the direction of the beam during the measurement 421 ] . 

The ion-optics of the Grand Raiden spectrometer is such that the particles scattered 
from the target position are focused vertically and horizontally at the focal plane. Using 
this property, the instrumental background was completely eliminated. While inelastically 
scattered a particles are focused at the focal plane, background events due to the rescattering 
of a particles from the wall and pole surfaces of the spectrometer show a flat distribution in 
the vertical-position spectra at the focal plane, as shown in Fig. [U The peak region in the 
vertical position spectrum was treated as true+background events. The off-center regions 
were treated as background only. Figure [2](a) shows the horizontal position spectrum for the 
112 Sn(a, a') reaction at 0°. The background spectrum has no distinct structure in the giant 
resonance region. Finally, we have obtained clean spectra by subtracting the background 
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spectrum from the true+background spectrum, as shown in Fig. 12(b). 




FIG. 1: Vertical-position spectrum at the focal plane of the Grand Raiden spectrometer, taken 
at 2.5°. The central densely-hatched region represents true+background events. The off-center 
sparsely-hatched regions represent only background events. The real events were obtained by 
subtracting background events from the true+background events. 



E 2500 -r 




Horizontal positionat focal plan (mm) 

FIG. 2: (a) Horizontal-position spectrum of the 112 Sn(a,o') reaction at 0°. The hatched region is 
background events, (b) Background-free spectrum 

The background-free "0°" inelastic spectra for the Sn isotopes are presented in Fig. [3j 
In all cases, the spectrum is dominated by the ISGMR+ISGQR peak near E x ~ 15 MeV. 
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There is an underlying continuum in the high excitation-energy region in the spectrum; it is 
reasonable to assume that this continuum, remaining after elimination of the instrumental 
background, is primarily due to contributions from excitation of the higher multipoles and 



quasifree knockout processes 



44|. 




FIG. 3: Excitation-energy spectra obtained from inelastic a scattering at 6i a b = 0.69° for all even-A 
Sn isotopes. 



III. DATA ANALYSIS 

We have employed the MDA procedure j^j] to extract the strengths of the isoscalar giant 
monopole resonance (ISGMR), the isoscalar giant dipole resonance (ISGDR), the isoscalar 
giant quadrupole resonance (ISGQR), and the high-energy octupole resonance (HEOR) in 
the Sn isotopes. The cross-section data were binned into 1-MeV energy intervals to reduce 
the statistical fluctuations. For each excitation energy bin from 8.5 MeV to 31.5 MeV, the 
experimental 17-point angular distribution (0 cm , E x ) has been fitted by means of the 
least-square method with the linear combination of calculated distributions -^-(9 cm , E x ) 
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defined by: 



da exp 7 do c £ l 

^ (9cm, Ex) — ^2<il{E x ) x — (6 cm ,E x ) (1) 



where -^-(9 cm , E x ) is the calculated distorted-wave Born approximation (DWBA) cross 
section corresponding to 100% energy-weighted sum rule (EWSR) for the L-th multipole. 
The fractions of the EWSR, cil{E x ), for various multipole components were determined by 
minimizing x 2 - This procedure is justified since the angular distributions are well charac- 
terized by the transferred angular momentum L, according to the DWBA calculations for 
a scattering. It was confirmed that the MDA fits were not affected by including L >7 



The DWBA calculations were performed following the method of Satchler and Khoa [451 ]. 
using the density-dependent single-folding model for the real part, obtained with a Gaussian 
a-nucleon potential, and a phenomenological Woods-Saxon potential for the imaginary term. 
Therefore, the a-nucleus interaction is given by: 

U(r) = V F (r) + iW/(l + exp((r - i2j)/fl/)) (2) 

where Vp(R) is the real single-folding potential obtained by folding the ground-state density 
with the density- dependent a-nucleon interaction: 

v D D G (r,r\p) = -v(l - /3p(r') 2/3 )exp(-|r - r'\ 2 /t 2 )) (3) 

where VuDG(f,r', p) is the density-dependent a-nucleon interaction, \r — r'\ is the distance 
between the center of mass of the a particle and a target nucleon, p(r') is the ground-state 
density of the target nucleus at the position r' of the target nucleon, /3=1.9 fm 2 , and i=1.88 
fm. W is the depth of the Woods-Saxon type imaginary part of the potential, with the 
reduced radius Ri and diffuseness aj. 

These calculations were performed with the computer code PTOLEMY Q]; with the 



input values modified 48) to take into account the correct relativistic kinematics. The shape 



of the real part of the potential and the form factor for PTOLEMY were obtained using 



the codes SDOLFIN and DOLFIN |49j. We use t 



re transition densities and sum rules for 



various multipolarities described in Refs. 13 



50 



51] . The radial moments were obtained 



by numerical integration of the Fermi mass distribution with the parameter values from 
Ref. [52J (listed in Table U). 



8 



TABLE I: Fermi-distribution parameters from Ref. 52]] . "c" is the adjusted half-density radius for 
the charge distribution and "a" is the diffuseness parameter. 



Target 112 Sn 114 Sn 116 Sn 118 Sn 120 Sn 122 Sn 124 Sn 



c(/m) 


5.3714 


5.3943 


5.4173 


5.4391 


5.4588 


5.4761 


5.4907 


a(/m) 


0.523 


0.523 


0.523 


0.523 


0.523 


0.523 


0.523 




FIG. 4: (a) Ratio of the elastic a-scattering cross sections to the Rutherford cross sections for 
112 Sn at 386 MeV. (b) Differential cross sections for excitation of the 2f state in 112 Sn. The solid 
lines are the results of the folding-model calculations. 




FIG. 5: Same as Fig. [U except for 120 Sn. 
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FIG. 6: Same as Fig. [H except for Sn 



TABLE II: OM parameters obtained by fitting elastic scattering data. Also listed are the B(E2) 



values for the corresponding 2f states from Refs. 



53 



54]. 



Target 


V (MeV) 


W (MeV) 


aj (fm) 


(fm) 


B(E2) (e 2 6 2 ) 


112 Sn 


33.9 


31.7 


0.60 


1.02 


0.24 


120 Sn 


33.4 


33.0 


0.63 


1.01 


0.20 


124 Sn 


34.0 


33.5 


0.61 


1.02 


0.17 



The optical-model (OM) parameters viz. the real part of the potential (Vp(r)), the 
Woods-Saxon type imaginary part of potential (W), the reduced radius (Ri), and the dif- 
fuseness (ay) in Eq. [2] were determined by fitting the differential cross sections of elastic a 
scattering measured for 112 Sn, 120 Sn, and Sn in a companion experiment; the results are 
listed in Table ITT1 The OM fits to the elastic scattering data for 112 Sn, 120 Sn, and 124 Sn, are 
shown in Figs. IHa),[5]^a), and[6](a), respectively. To test the efficacy of the OM parameters, 
DWBA calculations were carried out for the first 2 + states in these nuclei using a collective 
form factor and previously-established B(E2) values obtained from Refs. 53|, |54| (also listed 
in Table HT1) . Figures. IH(b), EJ^b), and[6]^b) compare the results of these calculations with the 
experimental data; indeed, the DWBA calculations reproduce the experimental differential 
cross sections for the 2f states well without any normalization. 

The contribution of the IVGDR excitation to the measured cross sections was subtracted 
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prior to multipole decomposition. Cross sections for exciting the IVGDR were obtained 
with DWBA calculations on the basis of the Goldhaber- Teller model and using the strength 
distribution obtained from photonuclear work |55j ]. 

Figs. [7] and [S] show the MDA fits to the experimental angular distributions of the differ- 
ential cross sections for the 16.5-MeV and 25.5-MeV energy bins in the inelastic-scattering 
spectra of 112 Sn and 124 Sn, respectively, along with the contributions from the L=0, 1 and 
2 multipoles. The ISGMR contribution is dominant in comparison to the other multipoles 
at E,j.=16.5 MeV. On the other hand, the ISGDR is the dominant contributor at E x =25.5 
MeV. 




FIG. 7: Angular distribution of 1-MeV bins centered at £^=16. 5 MeV for 112 Sn(a,a') and 
124 Sn(a,a'). The solid squares are the experimental data and the solid lines are the MDA fits 
to the data. Also shown are the contributions to the fits from L=0 (dashed line), L=l (dotted 
line), L=2 (dash-dotted line) and L=3 (small-dashed line) multipoles, as well as from the IVGDR 
(dash-dot-dotted line). 



IV. RESULTS AND DISCUSSION 



We have extracted strength distributions for L=0, 1, 2, and 3 multipoles over the energy 
range 8.5 MeV-31.5 MeV in all the Sn isotopes investigated in this work. These are displayed 
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FIG. 8: Same as Fig. [7] except for the 25.5-MeV energy bin (see text). 

in Figs. El [TUl dU and [I2l respectively. The strengths are related to the coefficients in 
Eq. 1 as (see, Refs. [13 



fill): 



9h A < r 2 > 

S (E X ) = 1 a (£.) (4) 



SlW = t^P4^(H < r 4 > -f < r 2 > 2 -106 < r 2 >)a 1 (E :E ) (5) 



5l> 2 (^) = ^4rL(2L + I) 2 < r 2L - 2 > a 2 {E x ) (6) 



where m, A and < r N > are the nucleon mass, the mass number, and the iVth moment of 

-2 



the ground-state density, and e=(4/ E2+B/ Eq)H /3mA; Eq and E2 are the centroid energies 
of the ISGMR and the ISGQR and have been taken as 80 A" 1 / 3 MeV and 64 A^ 1 / 3 MeV, 
respectively. 

It should be noted that although we employed calculated DWBA cross sections with up 
to L=7 in the MDA fitting procedure, it was not possible to reliably extract the strength 
distributions for L >4 because of the limited angular range (0°-9°). Further, there is a small, 
near-constant ISGMR and ISGQR strength up to the highest excitation energies measured 
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FIG. 9: (Color online) ISGMR strength distributions obtained for the Sn isotopes in the present 
experiment. Error bars represent the uncertainties from fitting the angular distributions in the 
MDA procedure. The solid lines show Lorentzian fits to the data. 

in this experiment. The raison d'etre of this extra strength is not quite well understood, 
iowever, similarly enhanced El strengths at high excitation energies were noted previously 



29 



30] and have been attributed to contributions to the continuum from three-body chan- 



nels, such as knockout reactions 



44l | . These processes are implicitly included in the MDA 



as background and may lead to spurious contributions to the extracted multipole strengths 
at higher energies where the associated cross sections are very small. This conjecture is sup- 
ported by measurements of proton decay from the ISGDR at backward angles wherein no 



56- 



such spurious strength is observed in spectra in coincidence with the decay protons [31 
581 ]; quasifree knockout results in protons that are forward peaked. A similar increase in the 
ISGMR strength at high excitation energies was reported as well by the TAMU group in 12 C 
when they carried out MDA without subtracting the continuum from the excitation-energy 



spectra 



18| 



The L = strength distributions were fitted with a Lorentzian function to determine 
the centroid energies and widths of the ISGMR. These fits are shown superimposed in 



13 







20 



10 



I 

30 



E (MeV) 



FIG. 10: (Color online) ISGDR strength distributions obtained for the Sn isotopes in the present 
experiment. Error bars represent the uncertainties from fitting the angular distributions in the 
MDA procedure. The solid lines show Lorentzian fits to the data. 

Fig. |9] the corresponding fitting parameters are presented in Table III II and compared with 



results from TAMU 
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211 ] . In this and subsequent comparisons and discussion, we refer 



only to the recent TAMU results because those are from comparable data and analysis — all 
other previous results on the Sn isotopes were from peak-fitting analyses of data taken at 
significantly lower energies. 

In order to compare with the available theoretical results, various moment ratios for the 
experimental ISGMR strength distributions have been calculated over the excitation-energy 
range, E x = 10.5-20.5 MeV, encompassing the ISGMR peak. The results are listed in 
Table IIVI The reasons for the difference between the present results and those from TAMU 
for 112 Sn and 124 Sn are not readily apparent but might be attributable to the fact that in 
their analysis the multipole decomposition is carried out after subtracting a "background" 
from the excitation-energy spectrum, whereas, as pointed out earlier, no such subtraction is 
required in the present analysis since the Sn(a, a') spectra obtained in our work have been 
rendered free of all instrumental background events. 
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FIG. 11: (Color online) ISGQR strength distributions obtained for the Sn isotopes in the present 
experiment. Error bars represent the uncertainties from fitting the angular distributions in the 
MDA procedure. The solid lines show Lorentzian fits to the data. 

Figure [10 shows the strength distributions of ISGDR. We observe a "bi-modal" distri- 
bution between ^=8.5 MeV and E x =31.5 MeV. This bi-modal pattern for the ISGDR 
has been observed in all nuclei investigated so far, both in the RCNP and TAMU mea- 
surements. This "low-energy" isoscalar L=l strength (LE) has engendered considerable 
interest and argument over the past few years. It is present in nearly all of the recent 
theoretical calculations in some form or the other, and at similar energies, although with 
varying strength. It has been shown 5JJ, |60] that the centroid of this component of the L=l 
strength is independent of the nuclear incompressibility and while the exact nature of this 



component is not fully understood yet, suggestions 
this component might represent the "toroidal" 60, 



lave been extended to the effect that 



61] or the "vortex" modes 62]. It is 



impossible to distinguish between the competing possibilities based on currently-available 
data [31] . There is general agreement, however, that only the high-energy (HE) component 
of this bi-modal distribution needs to be considered in obtaining a value of from the 
energy of the ISGDR. The strength distributions of the ISGDR, therefore, have been fitted 
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FIG. 12: (Color online) HEOR strength distributions obtained for the Sn isotopes in the present 
experiment. Error bars represent the uncertainties from fitting the angular distributions in the 
MDA procedure. The solid lines show Lorentzian fits to the data. 

with a two-Lorentzian function and the fitting parameters for the LE- and HE-components 
are presented in Tables [V] and |VI| respectively. It may be noted that because of the "spu- 
rious" strength at the higher excitation energies mentioned previously, the numbers for the 
extracted EWSR are significantly larger than 100% in some cases. 

The strength distributions of the ISGQR are shown in Fig. [Til These too were fitted with 
a Lorentzian function to determine the centroid energies and the widths. The fit parameters 
are presented in Table IVIIl 

The L=3 strength distributions (Fig. [T2|) show an enhanced strength at E x < 10 MeV. 
This part is, most likely, from the low-energy octupole resonance (LEOR). The LEOR 
represents the lhu component of the L=3 strength and has been reported in a number of 



nuclei previously [63|, |64|. The strength distributions were, therefore, fitted with a two- 
Lorentzian function to determine the centroid energy of HEOR (the high-excitation-energy 
component). The extracted HEOR peak-energies are presented in Table IVTlTI 
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TABLE III: Lorentzian-fit parameters for the ISGMR strength distributions in the Sn isotopes, as 
extracted from MDA. The quoted EWSR values are from the fitted Lorentzians. The results from 
TAMU work (from Gaussian fits), where available, are provided for comparison [l9. \^\. 



Target 


Eisgmr (MeV) 


r (MeV) 


EWSR a 


Reference 


112 Sn 
114 Sn 


16.1 ±0.1 
15-671^ 
15.9 ±0.1 


4.0 ±0.4 

e i o+0.40 
°- io -0.04 

4.1 ±0.4 


0.92 ± 0.04 
1 in+ a15 

i - iu -0.12 

1.04 ±0.06 


This work 
TAMU 
This work 


116 Sn 


15.8 ±0.1 


4.1 ±0.3 


0.99 ± 0.05 


This work 


118 Sn 


15.6 ±0.1 


4.3 ±0.4 


0.95 ±0.05 


This work 


120 Sn 


15.4 ±0.2 


4.9 ±0.5 


1.08 ±0.07 


This work 


122 Sn 


15.0 ±0.2 


4.4 ±0.4 


1.06 ±0.05 


This work 


124 Sn 


14.8 ±0.2 

■.K o 4 +0.13 


4.5 ±0.5 

5 00+ 013 
°- uu -0.53 


1.05 ±0.06 
1 06+ 010 

l.UO_ 20 


This work 
TAMU 


°Only statistical uncertainties are included; systematic errors, mostly from DWBA calculations, are ~15%. 



The primary focus of this work has been on the ISGMR because of its direct connection 
with the nuclear incompressibility. The excitation energy of the ISGMR is expressed in the 



scaling model [70[ as: 



Eisgmr = h\ % — (7) 

V m < r 1 > 

where m is the nucleon mass, < r 2 > the ground-state mean-square radius, and the 
incompressibility of the nucleus. 

The moment ratios, mi /mo, for the ISGMR strengths in the Sn isotopes are shown in 
Jig. [TBI and compared with recent theoretical results from Colo et al. (non-relativistic) 



65 



66( and Piekarewicz (relativistic) [67]. The calculations overestimate the experimental 
ISGMR energies significantly (by almost 1 MeV in case of the higher-A isotopes!). This 
difference is very surprising since the interactions used in these calculations are those that 
very closely reproduce the ISGMR energies in 208 Pb and 90 Zr. Indeed, this disagreement 



leaves open a puzzling question: Why are the tin isotopes so soft [67j? Are there any nuclear 
structure effects that need to be taken into account to describe the ISGMR energies in the 
Sn isotopes? Or, more provocatively, do the ISGMR energies depend on something more 
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TABLE IV: Various moment ratios for the ISGMR strength distributions in the Sn isotopes. All 
moments have been calculated over E x = 10.5-20.5 MeV. The quoted EWSR values are from the 
strength observed within this energy range. The results from TAMU work, where available, are 
provided for comparison 
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21] 



Target 


^ (MeV) 


( MeV ) 


( Me V) 


EWSR a 


Reference 


112 Sn 
114 Sn 


16.2 ±0.1 
16.1 ±0.1 


16.7 ±0.2 
16.05i°; 2 ^ 
16.5 ±0.2 


16.1 ±0.1 

■.k 90+0.10 
io.zo_ 10 

15.9 + 0.1 


0.73 ± 0.04 
0.86 + 0.05 


This work 
TAMU 
This work 


116 Sn 


15.8 ±0.1 
15.85 ±0.20 


16.3 ±0.2 


15.7 ±0.1 


0.86 ± 0.05 
1.12 ±0.15 


This work 
TAMU 


118 Sn 


15.8 ±0.1 


16.3 ±0.1 


15.6 ±0.1 


0.73 ± 0.04 


This work 


120 Sn 


15.7 ±0.1 


16.2 ±0.2 


15.5 ±0.1 


0.78 ± 0.05 


This work 


122 Sn 


15.4 ±0.1 


15.9 ±0.2 


15.2 ±0.1 


0.85 ±0.05 


This work 


124 Sn 


15.3 ±0.1 

14.50i°; 14 


15.8 ±0.1 
14.96l°i° 


15.1 ±0.1 

14 QQ+0.17 


0.77 ±0.05 


This work 
TAMU 



°Only statistical uncertainties are included; systematic errors, mostly from DWBA calculations, are ~15%. 



than the nuclear incompressibility, requiring a modification of the scaling relationship given 
in Eq. [_? 

There have been several attempts to explain this anomaly. One of the earliest was by 
Civitarese et al. to estimate the effect of pairing correlations on the energy of the 
ISGMR. The shifts obtained for the ISGMR energies of 100-150 keV across the Sn isotopic 
chain were insufficient to explain the experimental data. Piekarewicz and Centelles 72] have 
constructed a hybrid model having the same incompressibility coefficient [K oo =2?>0 MeV) 



as the FSUGold _3j while preserving the stiff symmetry energy of NL3 74J. This results 
in a considerably softer incompressibility coefficient for neutron-rich matter and produces 
a significant improvement in agreement with the experimental data on the ISGMR's in the 
Sn isotopes. However, as the authors point out, while the improvement in case of the Sn 
isotopes is unquestionable, an important problem remains: the hybrid model underestimates 
the ISGMR centroid energy in 208 Pb by almost 1 MeV, suggesting that the rapid softening 
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TABLE V: Lorentzian-fit parameters for the low-energy component of ISGDR strength distribu- 
tions in the Sn isotopes, as extracted from MDA. The results from TAMU work, where available, 
are provided for comparison Q] • 



Target 


Ele-isgdr (MeV) 


T (MeV) 


Reference 


112 Sn 
114 Sn 


15.4 ±0.1 
14.92±°;it 
15.0 ±0.1 


4.9 ±0.5 

o o 9 +0.26 

5.6 ±0.5 


This work 
TAMU 
This work 


116 Sn 


14.9 ±0.1 


5.9 ±0.5 


This work 




14.38 ±0.25 


5.84 ± 0.30 


TAMU 


118 Sn 


14.8 ±0.1 


6.1 ±0.3 


This work 


120 Sn 


14.7 ±0.1 


5.9 ±0.3 


This work 


122 Sn 


14.4 ±0.1 


6.7 ±0.3 


This work 


124 Sn 


14.3 ±0.1 

15 91+0.15 

io.oi_ 15 


6.6 ±0.3 
6-60±8;lI 


This work 
TAMU 



with neutron excess predicted by this hybrid model might be unrealistic. They also suggest 
that the failure of the FSUGold to reproduce the ISGMR energies might be due to missing 
physics unrelated to the incompressibility of neutron-rich nuclear matter; as an example of 
such missing physics, they mention the superfluid character of the Sn isotopes resulting from 
their open-shell structure. 

Calculations have also become available recently from the RMF approach with the DD- 
ME2 interaction 75], and these reproduce the centroids of the ISGMR in the Sn isotopes 



rather well 



68]. It is also seen that the DD-ME2 interaction falls within the constraints 



imposed by the experimental and K T values (see discussion below). Some concern 
has been expressed, however, that this agreement of the centroid energies might be just a 
coincidence since the ISGMR strength distributions for the Sn isotopes from this work a ppe ar 
to be not significantly different from those obtained from, for example, the FSUGold 76|. 

In calculations using the T5 Skyrme interaction within the quasiparticle time blocking 
approximation (QTBA) approach, Tselyaev et al. 69j have obtained the ISGMR strength 
distributions in all the Sn isotopes in good agreement with the experimental data, including 
the resonance widths. However, T5 has the associated value of only 202 MeV, which is 
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TABLE VI: Lorentzian-fit parameters for the high-energy component of ISGDR strength distribu- 
tions in the Sn isotopes, as extracted from MDA. The results from TAMU work, where available, 
are provided for comparison [lj]) Q] • 



Target 



Ehe-isgdr (MeV) 



r (MeV) 



EWSR 



Reference 



112 



Sn 



114 Sn 
116 Sn 

118 Sn 
120 Sn 
122 Sn 
124 Sn 



26.2 ±0.8 



26.28 



-0.32 
-0.23 



26.1 ±0.8 
25.9 ±0.6 
25.50 ± 0.60 
26.0 ±0.3 
26.0 ±0.4 
26.3 ±0.2 
25.7 ±0.5 



16.3 ±4.0 



10.82 



-0.39 
-0.36 



25.06 



+0.22 
-0.21 



13.9 ±3.4 
13.1 ±4.2 
12.00 ± 0.60 
13.1 ±2.0 

13.1 ± 1.9 
12.4 ± 1.1 

10.2 ± 1.6 
13-87181 



1.02 ±0.03 

7n+ 010 
u - ' u -o.io 

1.23 ±0.03 
1.02 ±0.03 
0.88 ±0.20 
1.20 ±0.03 
1.50 ±0.03 
1.47 ±0.03 
1.29 ±0.06 

o-93±8:i! 



This work 
TAMU 
This work 
This work 
TAMU 
This work 
This work 
This work 
This work 
TAMU 



"Only statistical uncertainties are included; systematic errors, mostly from DWBA calculations and the 
contributions at the highest energies (see text), are ~30%. 

significantly lower than that extracted earlier from the ISGMR's in 208 Pb and 90 Zr. While 
the agreement with the experimental data is impressive (and, indeed, reproduces the A- 
dependence rather well), it does leave the question of "softness" of the Sn nuclei unanswered. 
As the authors themselves state, the goal of their work has not been to solve the problem 
of the nuclear-matter incompressibility but to find under which conditions one can obtain 
reasonable description of the experimental data for the considered tin isotopes. 

The "superfluid" character of the Sn isotopes, resulting from pairing correlations in open- 
shell nuclei, has been investigated by Li et al. (77)]. In a self-consistent QRPA model 
that employs the canonical HFB basis and an energy- density functional with a Skyrme 
mean-field part and density-dependent pairing, they calculated the energy of the ISGMR 
for the Sn isotopes and looked at the effects of different kinds of pairing forces (volume, 
surface, and mixed). They find that, compared with the HF±RPA and HF-BCS-QRPA 
formalisms, the HFB±QRPA calculations lead to energies for the ISGMR in Sn isotopes 
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TABLE VII: Lorentzian-fit parameters of ISGQR strength distributions in the Sn isotopes, as 
extracted from MDA. The results from TAMU work, where available, are provided for compari- 
son 19, 3- 



Target 


Eisgqr (MeV) 


T (MeV) 


EWSR a 


Reference 


112 Sn 
114 Sn 


13.4 ±0.1 
13.48±g:it 
13.2 ±0.1 


7.0 ±0.5 

4 Q0+ 22 
4.yu_ 27 

6.8 ±0.4 


1.08 ±0.04 
0.88^ 
1.25 ±0.05 


This work 
TAMU 
This work 


116 Sn 


13.1 ±0.1 


6.4 ±0.4 


1.12 ±0.04 


This work 




13.50 ±0.35 


5.00 ±0.30 


1.08 ±0.12 


TAMU 


118 Sn 


13.1 ±0.1 


6.6 ±0.3 


1.08 ±0.03 


This work 


120 Sn 


12.9 ±0.1 


7.0 ±0.7 


1.04 ±0.04 


This work 


122 Sn 


12.8 ±0.1 


7.8 ±0.6 


1.25 ±0.04 


This work 


124 Sn 


12.6 ±0.1 

12 72+ - 11 
iZ - ' z -0.11 


7.7 ±0.9 

A on+0.32 


1.13 ±0.04 

0.8918:8 


This work 
TAMU 



"Only statistical uncertainties are included; systematic errors, mostly from DWBA calculations, are ^20%. 



that are significantly closer to the experimental values, in particular with the surface pairing 
forces and the SKM* interaction (Koo ~ 215 MeV) 78J]. Thus, while pairing effects lower the 
ISGMR excitation energies, one still needs to reduce the value by ~10% for achieving 
a reasonable agreement with the experimental data. 

A very intriguing possibility in explaining the "softness" of the Sn isotopes has been 



offered very recently by Khan 79|, l80[. The author asserts that, in analogy with the 



mutually-enhanced-magicity (MEM) effect observed in predictions of masses with different 



energy-density functionals 



81 



82| . the ISGMR energy in the doubly-magic nuclei might be 



anomalously higher. The obvious implication is that the calculations using interactions that 
are successful in describing the ISGMR in the doubly-magic nucleus 208 Pb would necessarily 
overestimate the ISGMR energies in the open-shell nuclei. If this effect is manifested in any 
significant way, the energy of the ISGMR in the non-doubly-magic Pb isotopes, 204 Pb and 
206 Pb, would be measurably lower than that in 208 Pb |8o| . In the only measurement of the 
ISGMR in 206 Pb reported so far [83|, this conjecture does not appear to hold. Still, precise 
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TABLE VIII: Lorentzian-fit parameters of HEOR strength distributions in the Sn isotopes, as 
extracted from MDA. The results from TAMU work, where available, are provided for compari- 
son 19 



u JU L 

1 



Target 



Eheor (MeV) 



T (MeV) 



Reference 



112 



Sn 



22.7 ±0.7 



7.2 ± 1.9 



20.63 



-0.30 
-0.28 



3.21 



1-0.30 
-0.28 



114 Sn 
116 Sn 

118 Sn 
120 Sn 
122 Sn 
124 Sn 



22.7 ±0.7 
22.3 ±0.6 
23.3 ±0.8 
22.1 ±0.6 
22.3 ±0.6 
22.1 ±0.6 
22.1 ±0.5 
19.12183 



7.2 ±2.1 
7.6 ± 1.7 
10.9 ±0.6 
5.9 ± 1.5 

7.5 ± 1.8 

5.6 ± 1.5 
8.1 ± 1.5 



3.30; 



-0.17 
-0.05 



This work 
TAMU" 
This work 
This work 
TAMU 
This work 
This work 
This work 
This work 
TAMU a 



"(mi/mo) ratios. 

measurements of the ISGMR in the Pb isotopes, using background-free inelastic spectra 
with high-energy a beams, would be worthwhile to fully examine this possibility. 
The incompressibility of a nucleus, K A , may be expressed as: 



K A ~ K vol (l + cA- 1 '*) + K T ({N - Z)/A) 2 + K Coul Z 2 A~^ z 



(8) 



Here, c ~ -1[84J, and K Co ui is essentially model independent (in the sense that the deviations 
from one theoretical model to another are quite small), so that the associated term can be 
calculated for a given isotope. Thus, for a series of isotopes, the difference K A — K C ouiZ 2 A~^ 3 
may be approximated to have a quadratic relationship with the asymmetry parameter, of 
the type y = A + Bx 2 , with K T being the coefficient, B, of the quadratic term. It has been 
established previously 85, Is6| that direct fits to the Eq. 8 do not provide good constraints 



on the value of K^. However, this expression is being used here not to obtain a value for 
Kqo, but, rather, only to demonstrate the approximately quadratic relationship between 
and the asymmetry parameter. 

From such an analysis of the ISGMR data in the Sn isotopes, we have obtained a value of 
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FIG. 13: (Color online) Systematics of the moment ratios mi /mo for the ISGMR strength distribu- 
tions in the Sn isotopes. The experimental results (filled squares) are compared with results from 



non-relativistic RPA calculations (without pairing) by Colo et al. 



calculations of Piekarewicz 



65, 



66J (filled circles); relativistic 



67] (triangles) ; RMF calculations from Vretenar et al. 68] (diamonds) ; 



and, QTBA calculations from the Jiilich group 



and Sn reported by the TAMU group 
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691 ] (sideways triangles). Results for Sn, "Sn 



211 ] are also shown (inverted triangles). 



K T = —550 ± 100 MeV (see Fig. 4 in Ref. [34]). This number is consistent with the value of 
K T = —370 ±120 MeV obtained from an analysis of the isotopic transport ratios in medium- 



energy heavy- ion reactions [87] . Incidentally, this value has 
—500 ±50 MeV that was quoted previously by this group 



Deen modified from the value of 
and referred to in Ref. 



34). 



It transpires that they had identified the quantity that they had obtained, K asy , as being 
identical to K T , the quantity that has been obtained from the ISGMR measurements; the two 



differ by a higher-order term [87 
obtained by Centelles et al. |q ll ] : 



90j. More recently, a value of K T 



-500 



-125 
-100 



MeV has been 



from constraints put by neutron-skin data from anti-protonic 
atoms across the mass table; here again, it would appear that what the authors have termed 
K T is actually the aforementioned K asy . Further, a value of K T = —500 ± 50 MeV has been 
obtained also by Sagawa et al. by comparing our Sn ISGMR data with calculations using 
different Skyrme Hamiltonians and RMF Lagrangians [92|. The K T value obtained from our 
ISGMR measurements has, thus, been verified by a number of different procedures involving 
quite different data. A more precise determination of K r will likely result from extending 
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the ISGMR measurements to longer isotopic chains. This provides strong motivation for 
measuring the ISGMR strength in unstable nuclei, a focus of current investigations at the 
new rare isotope beam facilities at RIKEN, GANIL, GSI, and NSCL 93-9j]]. 

Combined with the value of = 240 ± 10 MeV obtained from the ISGMR and IS- 
GDR data 
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65 



73 



971 ] , we now have "experimental" values of both and K r which, 



together, can provide a means of selecting the most appropriate of the interactions used 
in EOS calculations. For example, this combination of "experimental" values for Koo and 
K r essentially rules out a vast majority of the Skyrme-type interactions currently in use in 



nuclear stru cture calculations 
in Refs. [9qI. llQO] . 



90, 



98(. Similar conclusions were reached for EOS equations 



V. SUMMARY 



We have measured the strength distributions of the isoscalar giant resonances (ISGMR, 
ISGDR, ISGQR, and HEOR) in the even-A 112 ~ 124 Sn isotopes via inelastic scattering of 
386-MeV a particles at extremely forward angles, including 0°. The extracted parameters 
for these resonances are in good agreement with previously-obtained values where available. 
The ISGMR centroid energies are significantly lower than those predicted for these isotopes 
by recent calculations and point to the need for further theoretical exploration of applicable 
nuclear structure effects, especially the role of pairing in ISGMR strength calculations in the 
open-shell nuclei. The asymmetry-term, K T , in the expression for the nuclear incompress- 
ibility has been determined to be —550 ± 100 MeV from the ISGMR data in Sn isotopes and 
is found to be consistent with a number of indirectly extracted values for this parameter. 
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